Genuine Multipartite Entanglement Trends in Gapless-gapped Transitions of 

Quantum Spin Systems 
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We investigate the behavior of a genuine multiparticle entanglement, viz. the generalized geo- 
metric measure, in gapless-to-gapped quantum transitions of one-dimensional and two-dimensional 
spin models. The investigations are performed in the exactly solvable one-dimensional XY models, 
as well as two-dimensional frustrated Ji — J2 models, including the Shastry-Sutherland model. In 
particular, we find that multiparty entanglement acts as a better tool than the bipartite one in 
detecting such transitions in the frustrated quantum systems. 



I. INTRODUCTION 

Recent developments at the interface of quantum in- 
formation science and many-body physics indicate that 
quantum correlations, in particular, quantum entangle- 
ment [1] can potentially be a universal physical charac- 
teristic to investigate many-body phenomena. For ex- 
ample, entanglement has been proposed as a detector of 
quantum phase transitions in spin systems [2 4] and as 
a tool to develop and simplify efficient numerical simu- 
lations like density matrix rcnormalization group [5] . Of 
late, the developments involving the numerical simula- 
tions using matrix product states [6], projected entan- 
gled pair states [7] and tensor network states [8], have 
established a strong connection between entanglement 
and many-body theory. At the same time, the behav- 
ior of entanglement has been investigated in many-body 
systems like cold atomic gases in optical lattices [9], and 
trapped gaseous Bose- Einstein condensates [10]. Char- 
acterization of entanglement in such systems is mostly 
restricted to bipartite entanglement due to the general 
unavailability of computatable measures in the multi- 
partite scenario. However, in some cases, it turns out 
that the bipartite entanglement measures can not cap- 
ture the co-operative phenomena in the system [11, 12], 
and therefore, it is natural to look out for multipartite 
entanglement measures to investigate such many-body 
systems [13]. 

Multipartite entanglement measures, e.g. the geomet- 
ric measure [14] (cf. [15]) and global entanglement [16] 
have been used to describe many-body phenomena [17]. 
However, in general, they are hard to compute and it 
is therefore not possible to use them for states of arbi- 
trary many-body systems. Recently, a genuine multipar- 
ticle entanglement measure called generalized geometric 
measure (GGM) [18] has been introduced which can be 
easily computed for pure states in arbitrary dimensions 
and of an arbitrary number of particles. It has since been 
found to be useful to study the genuine multiparty entan- 
glement present in systems like resonating valence bond 
states and states of disordered systems [19]. In this pa- 
per, we apply the GGM to study quantum spin models, 
including frustrated ones [20, 21]. Frustration appears in 
a many-body system if it is not possible to simultaneously 
and independently minimize all the interaction terms of 



the corresponding Hamiltonian [21]. The characteriza- 
tion of such systems is typically hard to achieve [21]. We 
consider the following four classes of quantum spin sys- 
tems - 

1. the quantum one-dimensional (ID) XY model [22], 

2. the quasi-2D antifcrromagnctic Ji — J2 Hcisenberg 
model [23-25], 



3. the antifcrromagnctic J\ 
tice [20], and 



Jo model on a 2D lat- 



4. the Shastry-Sutherland model [26]. 

The choice of above models is due to their immense im- 
portance in understanding the different exotic phases in 
many-body systems including high-T c superconductiv- 
ity [27]. Moreover, the recent experimental realizations of 
such spin models in the laboratories [3, 28-30], for exam- 
ple, in optical lattice [3, 31], trapped ions [32], photons 
[33], and nuclear magnetic resonance [34], have led to 
the interesting possibility of the observation of the many- 
body effects described here in the laboratories. Towards 
the unfolding of such many-body effects, we apply the 
multiparty entanglement measure, the GGM, to detect 
the transition from gapless to gapped phases and vice- 
versa in these models. 

The paper is organized as follows: In Sec. II, we 
present a formal definition of the genuine multiparty en- 
tanglement measure, viz., the GGM, and we express the 
same in terms of easily computable Schmidt coefficients. 
The results are presented in Sec. Ill, where in each sub- 
section, we begin with a brief description of one of the 
quantum spin models. Sections III A and III B deal, re- 
spectively, with the anisotropic XY and the J\ — J2 mod- 
els, both in one dimension. In Sections IIIC and HID, 
we study two dimensional models - the 2D J\ — J2 model 
and the Shastry-Sutherland model. Finally, we draw our 
conclusions in Sec. IV. 



II. GENERALIZED GEOMETRIC MEASURE 

In this section, we present a brief description of the 
GGM, and show that it is efficiently computable for pure 
quantum states of an arbitrary number of parties. We 
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now define the genuine multipartite entanglement mea- 
sure, viz. the GGM. A pure quantum state AiA 2 ...A N , 
shared between N parties, A\,A%, ... , An, is said to be 
genuinely iV-party entangled, if it is not a product across 
any bipartite partition. The GGM of a pure quantum 
state \i/j)a 1 a 2 ...a n is defined as 

£{\'*1>)a 1 a*...a n ) = 1 - A^(\^)a 1 a 2 ...a n )- (1) 

Here A max (\ip) AlA2 .„ AN ) = ma,x\(<f>\ip} Al A 2 ...A N \ where 
the maximization is taken over all pure states \<f>) which 
arc not genuinely iV-party entangled. 

Let us enumerate some of the properties of GGM. 

1 . £ is non- vanishing for all genuine multiparty entan- 
gled states, and vanishing for others. 

2. £ is monotonically non-increasing under local 
(quantum) operations and classical communica- 
tion. (The proof follows from the theorem stated 
below and Ref. [35].) 

We now state a proposition, which shows how the GGM 
can be expressed in terms of Schmidt coefficients. 
Theorem: The generalized geometric measure (GGM) 
can be expressed as 

£(|V>)) = l-max.{\\. B \AUB = {A U ..., A N }, AT)B = 0}, 

(2) 

where Xa-.b is the maximal Schmidt coefficient of 
\4'}a 1 a 2 ...a n in the A : B bipartite split. 



III. GGM DETECTS GAPLESS - GAPPED 
QUANTUM TRANSITIONS 

In this section, we consider a series of paradigmatic 
quantum spin systems. They are taken up one by one in 
the different subsections. Each subsection begins with a 
brief description of the Hamiltonian corresponding to the 
quantum system under study. Subsequently, we study 
the behavior of the ground state of these models, and in- 
vestigate the advantage of considering the genuine mul- 
tipartite entanglement measure in these models. 



A. Anisotropic XY model 

The Hamiltonian for the one-dimensional anisotropic 
quantum XY model of N quantum spin-half particles, 
arranged in an ID array, is given by [22] 



Hxy = 



J 



f N 

£(l + 7)«+i 
\i=i 



N 



+ (l-7)o-W+i +fcE°* 



(3) 

where J is the coupling constant for the nearest neighbor 
interaction, 7 £ (0, 1] is the anisotropy parameter, cr's 
are the Pauli spin matrices, and h represents the exter- 
nal transverse magnetic field applied across the system. 



In all the models considered in this paper, we impose the 
periodic boundary condition. The quantum XY Hamil- 
tonian can be diagonalized by applying Jordan- Wigner, 
Fourier, and Bogoliubov transformations [22]. At zero 
temperature, it undergoes a quantum phase transition 
driven by the transverse magnetic field at A = h/J = 1. 
Moreover, it is also known that the model is gapped for 
all field strengths except at the point where the quan- 
tum phase transition occurs. Such transitions have been 
detected by using bipartite entanglement measures like 
concurrence and multipartite entanglement measure like 
geometric measure [17, 36]. We investigate the behavior 
of the genuine multipartite entanglement measure viz., 
the GGM, of the ground state, when it crosses from one 
gapped phase to another, through the gapless point. 

The ground state of the system represented by the 
quantum XY Hamiltonian, as given in Eq. (3), can be 
analytically obtained by using Majorana fermions and it 
is also possible to get the eigenvalues of the local states 
corresponding to the ground state, in different biparti- 
tions [4, 22] which are given by 



n 



1 +(-!)*<!/< 



0, 1 Vi, (4) 



where u^s are the eigenvalues of Gl- The latter is in turn 
given by 



G L = 



90 ■ ■ gL-i 
-9l-i ■ ■ go 



with 



91 = Y 



2 71 



cos 4> — A — 17 sin <; 
I cos — A — sin < 



(5) 



(6) 



In Fig. 1, we have plotted the GGM and the derivative 
of the GGM for the ground state of the XY model with 
respect to the transverse magnetic field, A, for different 
values of the anisotropic constant 7. The logarithmic di- 
vergence of the derivative of GGM captures the presence 
of the quantum phase transition at A = 1 . When 7=1, 
which corresponds to the Ising model, the genuine multi- 
partite entanglement is maximum when compared to the 
systems with lower values of 7. 



B. Quasi-2D Frustrated Ji J 2 Model 

We will now consider the frustrated quasi two- 
dimensional Ji — J 2 Hciscnbcrg model, in which the near- 
est neighbor couplings, J\, and the next-nearest neighbor 
couplings, J2, are both antiferromagnetic. It was found 
that solid state systems like SrCu02 can be described by 
this model [37]. Moreover, advances in the field of cold 
atomic systems promise to create and control such mod- 
els in the laboratory [3]. The Hamiltonian of this model, 
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FIG. 1. (Color online.) GGM of the transverse XY model. 
The GGM (actually £ — 1/2) and its derivative (both dimen- 
sionless) are plotted on the vertical axis for the anisotropic 
transverse XY model for different anisotropy parameters 7, 
against the dimensionless system parameter A on the hori- 
zontal axis. The dashed (blue), circled (pink), and dotted 
(green) lines are respectively for the Ising (7 = 1), 7 = 0.8, 
and 7 = 0.2 models. The derivatives of GGM diverges at the 
quantum critical point A = 1. 



with N lattice sites on a chain, can be written as 



N 



(7) 



i=i 



where J\ and J2 are both positive. In the parameter 
space, a = J2/J1 = 0.5 is known as the Majumdar- 
Ghosh point, and the system is highly frustrated there. 
For an even number of sites, the ground state at the 
Majumdar-Ghosh point is doubly degenerate, and the 
ground state manifold is spanned by the two dimcrs 

\tPmg) = n^ = / 1 2 (|0) 2 i]l>2t±i - |l)2i|0) 2i± i). Note that 
the model is gapped at this point [23]. For a = 0, 
the Hamiltonian reduces to the spin— g Hcisenberg an- 
tiferromagnet and hence the ground state, which is a 
spin fluid state having gapless excitations [38], can be 
studied by Bethe ansatz. At other points, the ground 
state and the energy gap of this model were considered 
by using exact diagonalization, density matrix renormal- 
ization group method, bosonization technique, etc [25]. 
It is known that at a w 0.2411, a phase transition from 
fluid to dimerization occurs [39] . In the weakly frustrated 
regime, < a < 0.24, the system is gapless, and therefore 
critical [23, 24]. The system enters a dimerized regime, 
and is gapped, for higher values of the coupling parame- 
ter. It is known that bipartite entanglement cannot de- 
tect this gapless to gapped transition [12]. 

The multisite entanglement measure, GGM, is how- 
ever able to detect the same. We perform exact diag- 
onalization of the Hamiltonian using TITPACK ver. 2 
developed by H. Nishimori [40] to find the ground state 
and then compute the GGM. In Fig. 2(a), a discontinuity 
in £ at a = 0.5 is visible, which signals the Majumdar- 



Ghosh point. It is also clear from Fig. 2(a) that the GGM 
changes its curvature at the gapless to gapped transition. 
To identify the exact phase transition point, we introduce 
an order parameter by using GGM (£), which we call the 
bound GGM, defined as 



£ 



d 2 £ 
da 2 



(8) 



We find that the fluid-dimer transition can be detected 
by the vanishing of the bound GGM at a w 0.23. The 
discontinuities in the GGM curves, apart from the dis- 
continuity at a = 0.5, are arguably due to the finite and 
small system sizes. Note that the GGM curves asymp- 
totically go to zero since the spin chain decouples into 
two spin chains with nearest neighbor interaction cou- 
plings, J2, for very high values of the driving parameter 
a. It may be noted here that although £b = is a sim- 
ple differential equation, its solution can be found only 
by knowing the relevant boundary conditions, which, in 
this case, can be thought of as being provided by the 
Hamiltonian corresponding to the system. This renders 
the differential equation intractable (analytically or nu- 
merically), and we do not follow that path. 
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FIG. 2. (Color online.) GGM and bound GGM for the 
quasi-2D frustrated antiferromagnet. In (a), the GGM (di- 
mensionless) is plotted on the vertical axis and system param- 
eter a (dimensionless) is plotted along horizontal axis. The 
change in curvature between a = 0.2 and a = 0.4 indicates 
the gapless to gapped phase transition. In (b), the bound 
GGM (dimensionless) is plotted on the vertical axis and sys- 
tem parameter a (dimensionless) is plotted along horizontal 
axis. The fluid-dimer transition is signaled by the vanishing 
of the bound GGM. 
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FIG. 3. (Color online.) GGM for the 2D frustrated antifer- 
romagnet. The GGM (vertical) is plotted against the system 
parameter a (horizontal). Both the quantities are dimension- 
less. 



C. 2D Frustrated Ji - J 2 Model 



obtain the ground states, we use exact diagonalization 
technique as mentioned in the preceding subsection. In 
Fig 3, we plot the GGM as a function of the driving pa- 
rameter a. From the behavior of the GGM, we predict 
that the dimer to collinear transition point is in the range 
a £ (0.61,0.62) for N = 16, and a ps 0.81 for N = 24. 
The indicator here is a maximum appearing between a 
long-span increase and a long-span decrease in the GGM 
curve with respect to a. Note that there are other short- 
span increases and decreases (with local maxima in be- 
tween), which are attributable to finite size errors. The 
Neel to dimer transition is estimated to occur at a ~ 0.3 
for both the cases, indicated by a sharp bend there of the 
GGM curve. Just like in the ID case, the GGM curves 
asymptotically go to zero for high a. This is because the 
spin lattice decouples into two spin lattices with nearest 
neighbor interaction couplings, J 2 , for very high values 
of the driving parameter a. 



We now consider an arrangement of quantum spin— i 
particles on a 2D square lattice, where nearest neighbor 
spins (both vertical and horizontal) on the lattice are cou- 
pled by Heisenberg interactions, with coupling strengths 
Ji , and where all diagonal spins are coupled by the same 
interactions, with coupling strengths J 2 - Both J\ and J 2 
are considered to be positive. The model has attracted 
a lot of attention [41] due to its connection with high 
T c -superconductors and its similarity with magnetic ma- 
terials like Li 2 VOSi0 4 and Li 2 VOGe0 4 [42]. Although 
the different phases of the ground state of this model have 
been predicted by different numerical as well as approxi- 
mate analytical methods [43] , some debates remain. The 
Hamiltonian of the system is given by 

H 2 d = Ji (Ji-Vj + Jz Vi-Vj, (9) 

(NN) (diagonals) 

where J\ and J 2 are antiferromagnetic. 

In the classical limit, only a first-order phase transi- 
tion from Neel to collinear at a e J2/J1 = 0.5 is ex- 
hibited by this model. The nature of the phase diagram 
changes when quantum fluctuations are present, and in 
this case, the exact phase boundaries are not known. 
Based on exact diagonalization, series expansion meth- 
ods, field-theory methods [43], etc., one expects that 
there are two long range ordered (LRO) ground state 
phases in the system, that are separated by quantum 
paramagnetic phases without LRO. These investigations 
predict that the first transition from Neel to dimer oc- 
curs at a = a\ £ (0.35, 0.45) while dimer to collinear 
transition happens at a = a 2 £ (0.62,0.68). There are 
recent experimental observations and proposals of detect- 
ing these phases in the laboratory, and they demand a 
precise quantification of the low temperature phase dia- 
gram of this model. 

We investigate the behavior of genuine multipartite en- 
tanglement of the ground state and study its effective- 
ness to detect the transitions present in the system. To 



D. The Shastry-Sutherland Model 
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FIG. 4. The Shastry-Sutherland lattice. The solid lines repre- 
sent the nearest neighbor interactions with coupling strength 
Ji between the lattice sites and the ones joined by the 
dashed lines represent next neighbor interactions with cou- 
pling strength J2. 

In this section, we study the entanglement properties 
of systems where the interaction between particles can 
be modeled by the Shastry-Sutherland Hamiltonian [26]. 
In the insulators like SrCu2(B03)2, the low-energy spin 
excitations reside on the spin-half copper ions which lie 
in two dimensional layers decoupled from each other. 
The antiferromagnetic exchange couplings between the 
Cu ions is identical to the Shastry-Sutherland Hamilto- 
nian. The lattice with schematic interactions, for this 
model, is given in Fig. 4 and the Hamiltonian is given by 

Hss = Ji'^Vi-Vj + h ^ (10) 

NN NNN 

Here, Ji(> 0) corresponds to nearest neighbor interac- 
tion (indicated by solid lines in Fig. 4) and J 2 (> 0) cor- 
responds to specific next nearest neighbors (indicated by 
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broken lines in Fig. 4). It is known that a simple product 
of singlet pairs, on the diagonal links, is the ground state 
of Hss for sufficiently large a = J2/ J\- It has been pre- 
viously reported [44, 45] that the system undergoes two 
quantum phase transitions driven by the quantum fluctu- 
ations: one is from Neel to an intermediate phase and the 
other one is from that intermediate phase to dimer. The 
nature of the intermediate phase is not yet clearly under- 
stood [46]. In Fig. 5, we plot the GGM as a function of 
a for 16 particles on a square lattice interacting via the 
Shastry-Suthcrland Hamiltonian. There are clear signa- 
tures of these phase transitions in the figure at a « 1.05 
and at a w 1.53, as have been predicted by other meth- 
ods [44, 45]. Due to the computational constraints we are 
able to report our finding only for N = 16. More specif- 
ically, we expect that studying the Shastry-Sutherland 
model requires an exact square structure with even num- 
bers of spins on each side, and hence the next relevant 
lattice size is N = 36. For N = 16, we have again used 
TITPACK ver.2 [40]. 
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FIG. 5. (Color online.) GGM for the Shastry-Sutherland an- 
tiferromagnet. GGM (on vertical axis) is plotted with respect 
to the system parameter a (on horizontal axis) for 16 spins. 
Both the quantities are dimensionless. 



IV. CONCLUSIONS 

An important classification scheme for multipartite en- 
tangled quantum states is according to their separability 
in different partitions. The complexity of such a classi- 
fication makes it difficult to obtain a unique multiparty 
entanglement measure, even for pure quantum states. A 
comparison with the situation for mixed bipartite states 
is relevant here. While the entanglement of pure bipar- 
tite states can be uniquely characterized by a single en- 
tanglement measure, a variety of different measures ex- 
ist for mixed bipartite states. In the case of multiparty 
quantum states, one can analogously have "pure" and 
"non-pure" forms of entanglement, corresponding respec- 
tively to genuine multiparty entanglement and ones that 
are not genuine. In this sense, the generalized geomet- 
ric measure, which is a measure of genuine multiparty 
entanglement, quantifies the pure form of multiparty en- 
tanglement that is present in a multiparty quantum state. 

We employed this measure to predict phase diagrams 
in quantum many-body systems. We began by using the 
measure to detect the quantum fluctuation driven phase 
transition in an exactly solvable model, viz. the quantum 
XY model. We subsequently applied the generalized ge- 
ometric measure to prototype frustrated quantum spin 
models, in the quasi two dimensional antiferromagnetic 
J\ — J2 model, the two dimensional antiferromagnetic 
J\ — J2 model and the Shastry-Sutherland model. The 
ground states and the corresponding phase diagrams, for 
the frustrated models, are not known exactly, although 
there have been several predictions by different meth- 
ods. We use exact diagonalization techniques to investi- 
gate the multipartite entanglement of the ground states 
of the frustrated models. We find that the genuine multi- 
partite entanglement measure is capable of distinguishing 
between the gapless and gapped phases in these models. 
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